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We present a combined analysis of SESAM and TxL data for the flavor singlet axial coupling G\ of the proton, 
which is very helpful to stabilize the disconnected signals at small quark masses. From connected and disconnected 
contributions we use the tadpole improved renormalization constant Za and obtain G\ = 0.21(12). 
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1. Introduction 

The flavor singlet axial coupling of the proton 
is defined as 

Sm'^a = {P\uit,lbU + ht^lbd + s7m75s|P), (1) 

where s^ denotes the components of the proton 
polarization vector. In the naive parton model 
G\ is related to the fraction of the proton spin 
carried by the quarks. From the measurement of 
the first moment of the spin dependent proton 
structure function, g^, in deep inelastic polarized 
muon proton scattering the EMC experiment pi 
found a small value. 



Gil -0.12(17) (2) 

which led to the "proton spin crisis" . The re- 
sult for G\ indicates that the contribution to the 
proton spin from the quarks is small. New exper- 
iments, including proton, neutron and deuteron 
data Ql, have shifted the value up to G\ — 
0.29(6) which is still far away from the Ellis- Jaffe 
QCD sum-rule prediction 0]: 
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0.579(25) 



(3) 



Albeit a lot of theoretical and experimental work 
has been done in the meantime, a clear under- 
standing of the proton spin is still lacking. In this 
context lattice calculations are very illuminating 
though hampered by large fluctuations Q. 



2. Lattice techniques 

On the lattice both connected and disconnected 
diagrams for G\ (see figure |l|) can be calculated 
from ratios of correlation functions M . 





(a) connected 



(b) disconnected 



Figure 1. Connected (a) and disconnected (b) di- 
agrams contribute to the flavor singlet axial cou- 
pling G\ of the proton (j^ = qji_,j5q)- 

The SESAM collaboration has analyzed 200 
16^ X 32 gauge configurations at each of their 
four quark masses corresponding to m.,:/mp = 
0.833. ..0.686. The configurations have been gen- 
erated previously in the standard Wilson dis- 
cretization scheme with Uf = 2 mass degener- 
ate quark fiavors, at (3 = 5.6. In addition 200 
TxL gauge configurations (24^ x 40 lattices) at 
the same coupling (3 = 5.6 have been analyzed. 
The quark mass of the TxL lattices corresponds 
to the lightest SESAM quark mass. Details of the 
simulations can be found in 
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2.1. Connected contributions 

For the amplitudes from the connected dia- 
grams (see figure |](a)) in the matrix element, eq. 
11 we have applied the global summation method 
and the insertion technique as described in H . 



From the ratio 

"EsiP^ J2y,yo [^TpTsg] (y, yo)P) 
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(X! [97^759] (y, 2/0)) , 



(4) 



where P is an interpolating operator for 
the proton, the connected amphtude Cq — 
{P\Qlfi'y5<l\P)conn Can be obtained from the 
asymptotic hnear slope 



i?Sr/M(^) t^oo 



A+{P\qj^j5q\P)t 



(5) 



2.2. Disconnected contributions 

For the disconnected amplitudes Dq — 
{P\qi tilbq\P) disc (see figure 0(b)) we have used 
the plateau accumulation method (PAM) H]: 
from the partial summation 
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Dq follows from the asymptotic time dependence 

{t —>■ 00) in 
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it,Ata,At) = B + Dq{t-At-Ata). (8) 



For the calculation of the axial vector quark 
loops we have used the spin explicit stochastic 
estimator technique with complex Z2 noise and 
100 stochastic estimates per spin component and 
gauge field configuration B. 



3. Raw data 



Figure g displays R^^ it) for the connected 
amplitudes Cu,d on the four SESAM quark 
masses. The results from the linear fits are sum- 
marized in table |l]. For the lightest SESAM quark 
mass {usea — 0.1575) we have also calculated Cu,d 
on the TxL lattice. Note that the volume effects 
are less than 7%. 
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Figure 2. R'^/^'^ 

configurations (spin-averaged) . The linear fits 

give the connected amplitudes C„ and Ca- 
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^U 


Cd 


0.1560 

0.1565 

0.1570 

0.1575(SESAM) 

0.1575(Txi) 


1.100(8) 

1.102(11) 

1.025(12) 

1.018(24) 

1.086(11) 


-0.307(3) 

-0.308(4) 

-0.297(8) 

-0.295(10) 

-0.310(6) 



Table 1 

Connected amplitude Cu,d from the linear fits to 
RT^''{i) (SESAM and TxL). 



In Figure I we plotted i?5^^'^(i) (At = Ato = 
1) for the disconnected amplitudes with Hgea = 
0.156 (SESAM) and n^ea = 0.1575 [TxL). Albeit 
the statistical errors for Dq are large we observe 
comparable quality of signals from both plots. 



4. Results 

The chiral extrapolations of Cu.d and Dqs to 
the light quark mass mught are shown in figure y. 
Since the finite volume effects are small we have 
taken the TxL data at Ksea = 0.1575 for the ex- 
trapolation in Dq. Within the errors the flavor 
symmetry Dq = Ds in the disconnected ampli- 
tudes is preserved and we use only the connected 
amplitudes for G\ and G\. 

With the Wilson discretization both the flavor 
singlet and non-singlet current need renormaliza- 
tion. We have used the first order (tadpole im- 
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Figure 4. Linear chiral extrapolations of Cu.d 
and Dq^g to the light quark mass mught- 



Figure 3. Ra^^' (t) for K^^a = 0.156 (SESAM) 
and Ksea — 0.1575 (TxL) (spin-averaged) . The 
linear fits give the disconnected amplitude Dq . 

proved) result from lattice perturbation theory 
for Z^ and Z^'^Ej. The renormalized coupling 
constants of the proton and the contributions Aq 
to the proton spin are listed in table 0. 
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0.62(7) 


-0.28(6) 


-0.12(7) 


G\ 


G\ 


G\ 


0.21(12) 


0.907(20) 


0.484(18) 



Table 2 
Renormalized coupling constants of the proton 
and the contributions Aq to the proton spin. 

5. Conclusion 

We have calculated connected and discon- 
nected contributions to the flavor singlet axial 
vector coupling of the proton in a full QCD 
nf — 2 lattice simulation with Wilson fermions. 
We have checked the volume effects for the con- 
nected amplitudes when we increase the lattice 



size at the same coupling (/3 = 5.6) and quark 
mass. In the analysis of the disconnected con- 
tributions we have replaced the SESAM results 
for the lightest quark mass by the TxL data and 
obtained a better signal for Dq. Nevertheless, 
the errors are dominated by the large statistical 
fluctuations in the disconnected amplitudes. We 
find G\ = 0.21(12) to be consistent with the re- 
sult from experiment and with previous quenched 
estimates ^^. For the triplet couphng we get 
G\ = 0.907(20) which is 30% smaller than the 
experimental value, G\ = 1.2670(35) g. To fix 
the systematic errors in the lattice calculation we 
would need a scaling analysis as well as a non- 
perturbative determination of Z^^ and Z^- q 
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